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Description 

Background of the Invention 

[0001 ] The invention relates to monitoring the operational performance of fluid storage systems 
[0002] Large quantities of liquids and similar materials are often stored in bulk storage containers or tanks which 
may be located above-ground, partially above-ground, or completely below ground. Such containers ortanks are gen- 
erally connected by piping to flow-meters or dispensers. 

[0003] For example, underground storage tanks (USTs) and. occasionally, above-ground storage tanks (AST's) are 
used to store petroleum products and fuel to be dispensed at automobile service stations, trucking terminals, automobile 
rental outlets, and sinriHar operations through gasoline, diesel, or kerosene dispensing pumps. Fuel product is generally 
dehyered to such facilities by a gravity drop from a compartment in a wheeled transport means such as a fuel delivery 
truck. ASrs or USTs are often located at central distribution locations so that product can be subsequently withdrawn 
from the tank system to be transported for delivery to a variety of such facilities. A distribution location with USTs or 
ASTs may receive deliveries of product from, e.g.. a pipeline spur, wheeled transport, a barge, or a rail car. 
[0004] Direct observation of the operating condition of such tanks and storage containers is difficult or impossible 
The vanous methods for identifying the amount of product in tank systems have varying levels of accuracy repeatability" 
and perfomiance. Moreover, the accuracy of devices which measure the amount of product dispensed from the storage 
containers and tanks differs greatly, and may or may not be temperature compensated. The amount of product actuallv 
de ivered to the tank system is often measured inaccurately and, frequently, not at all. Rather, the owner of operator 
of the tank or vessel usually records the invoiced amount of product delivered as the actual amount introduced to the 
tank system, without having any means of confimrilng whether the Invoiced amount of product delivered is correct 
[0005] Consequently, effective management of such facilities is complicated by the numerous errors in the various 
measuring devices and procedures used to establish a baseline for management, planning and decislonmakinq Ef- 
fective management requires the following: 

1 . Accurate measurement of the volume stored in the system. 

2. Accurate detemiination of the volume dispensed from the system. 

3. Accurate detemiination of the amount of product introduced Into the system. 

4. Identification of volumes added to or removed from the tank system which are not otherwise recorded 

5. Rapid identification of leakage from the tank system. 

6^ Continuous monitoring and diagnosis of the operating performance of all of the component measuring devices 
of the system. ^ 

7. Continuous analysis of sales data to predict demands of product from the system 

B. Detemiination of optimal reorder times and quantities as a function of ordering, transportation, holding and 
penalty costs in order to minimize total costs of operation and/or to maximize profits. 

[0006] Traditionally, these functions were perfomned crudely or, in many cases, not at all. Volume measurements 
were, and m many instances still are. based on imperfect knowledge of the geometry, dimensions, and configuration 
of the storage vessel. Also, dispensing meters are frequently miscalibrated. This is true even when tank systems are 
regulated, due to the breadth of tolerance permitted for individual sales as related to total tank volume For examole 
deliveries from the delivery vehicle are almost always unmetered, additions of product from defueling vehicles are 
typically undocumented, and theft of the product is not uncommon. 

[0007] Leakage of product has, in recent years, assumed a dimension far in excess of the mere loss of the product 
Environmenta daniage can, and frequently does, expose the operatorto very large liabilities from third party l^igation 
m addition to U.S Environmental Protection Agency (EPA).mandated remediation which can cost in the range of hun- 
dreds of thousands of dollars. The ERA'S requirements for leak detection are set forth in EPA Pub. No. 510-K-95-003 
Straight Talk On Tanks: Leak Detection Methods For P etroleum Underoround Storage Tanks and PIpino (July 1 991)' 

which IS incorporated herein by reference. ^—^ '' 

[0008] To address these concerns, Statistical Inventory Reconciliation (SIR) was developed. The SIR method con- 
sists of a computer-based procedure which identifies all of the sources of error noted above by statistical analysis of 
the vanous and unique pattems that are introduced into the Inventory data and. in particular, into the cumulative vari- 
ances in the data when viewed as functions of product height, sales volumes, and time. 

Summary of the Invention 

[0009] The present invention relates to an automatic SIR system that may continuously and automatically collect 
data from completely above-ground, partially above-ground, and completely below ground containers for statistical 
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analysis. The invention addresses a variety of physical, business, operational and envlronnnental issues associated 
with the bulk storage of liquids or pourable solids. 

[0010] The present invention is an application of SIR that greatly enhances the ability to manage a facility effectively. 
It provides the means to characterize exactly the geonnetry, dimensions, and configuration of the storage vessel, identify 
overages and shortages In deliveries and unexplained additions and removals of product, and provide an accurate 
assessment of overall dispensing meter calibration. In addition, by accounting for such discrepancies, the present 
invention pennits identification of leakage at rates less than .1 gallon per hour in all of its estimates to any prescribed 
tolerance; By increasing the number of measurements taken, the estimates can be derived at any desired level of 
tolerance. 

[0011] The method of the present invention makes no assumptions as to the precision of any of the measuring 
devices used in various system configurations. Precision and calibration accuracies are derived from the data alone. 
Also, it is not assumed that the tank system is leak free; the leak status of the system is determined from the data alone. 
[0012] The method derives tank geometry, dimensions, and configuration, and their impact on the totality of cumu- 
lative inventory variances, as a function of product height in the tank. Correctness of dispensing meter calibration is 
verified in a similar manner by testing for randomness of cumulative variances as a function of varying sales volumes. 
Having confimned that such remaining residual variances are random, reflecting only the inherent random noise of the 
measurement devices, the present method analyzes departures of the cumulative variance from the bounds determined 
by the calculated random noise level. All calculations as to the volumes added, removed, metered or leaking are based 
upon extended successive, simultaneous observations of meter and gauge readings. The number of observations 
incorporated in each such calculation is determined by computing confidence bands forthe parameters of interest and 
extending data collection as necessary to achieve predetermined tolerances. 

[0013] For example, the method of the present invention is capable of distinguishing between continuous losses 
consistent with leakage and one-time unexplained removals of the fluid product from the tank. The method may be 
used to ensure the accuracy of computed delivery volumes, which are detemnined and reported with confidence bound- 
aries calculated for estimated delivered quantities. 

[0014] The method can also be used to control and monitor the accuracy of purchase costs of fluids such as petroleum 
which are delivered to tanks. For example, motor fuel retailers may be charged by wholesalers for either net or gross 
volumes purported to have been delivered. A determination that purchase charges are appropriate thus requires fre- 
quent simultaneous readings of sales, tank volumes and temperatures, which can be accomplished using the method 
of the present invention. 

[0015] To accomplish these goals, the present invention involves estimating changes in product volume in a tank 
based on multiple data points and their respective likely errors measured continuously over a period of time, A software 
program is used to implement an algorithm that employs concepts from matrix theory and mathematical statistics. The 
algorithm includes generating the product of a matrix and its transpose by successive additions of partial products of 
partitions of the matrix and their corresponding transposed matrix partitions to minimize the storage requirements of 
the data collected. The compressed matrix data constitutes a complete and sufficient statistic for the parameters of 
interest. The algorithm thus permits the accumulation and storage of a large amount of data in a condensed form 
without sacrificing statistically useful information, to obtain a statistically significant result with the required, accuracy 
and reliability. 

[0016] Thus, one object of the present invention is to determine the accuracy and consistency of devices used to 
measure volume of product added to, removed from, and present in a fluid storage system. 

[0017] Another object of the invention is to identify and quantify additions of material to the system, but not recorded 
as such, and volumes of product removed from the system which are not registered by measuring devices or othenwise 
recorded. 

[0018] Another object of the invention is to discriminate between discrete one-time unrecorded removals of product 
from the system and continuous losses consistent with leakage. 

[0019] Another object of the invention is to identify and provide early warning of product leakage from all parts of the 
system, extending from the fill point to the point of discharge, and to confimn the validity of resulting leakage warnings, 
[0020] Another object of the invention is to detemnine secular, seasonal trends and repetitive special demands to 
provide short and long term estimates for demand of the product, and to provide optimal reorder quantities and a 
delivery schedule for the system. 

[0021] Afurther object of this invention is to accomplish alt of the foregoing in a fully automated system that requires 
no human intervention, other than as an option available to the operator to enter quantities of material Teportediy 
delivered for comparison with those computed. 

[0022] In general, in one aspect, the invention features a method of monitoring a fluid storage and dispensing system. 
The system has measurement apparatus for measuring a volume associated with the system. A plurality of measure- 
ment data is collected from the measurement apparatus in a form readable by a computer and stored in a compressed 
matrix format In a computer memory. The compressed matrix format is statistically analyzed to detemnine operational 
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monitoring information. 

n-lng «.e measurement e^L,.. under m^^^^TJll^r ""^ 

SerrT o, ..etemperatu're o„.e v^rrsu^ rrprxroTrnr.srrz^ = 

of maasuremant data Is collectad froiX m?ru™m.^t ? ™ associated with the systan,. A plurality 

cornprassad mathx (omiat In a con^X ImoTZ comTS ™, Jr "f"^".''^ ' i" » 

wllether a quantity ol (luid ,»mo»ad fi^n, SysTam la Zi^hvf 71 " •"^'V"". » deten„lne 

dataotad In the siitan, Is delivered ' °*'' '" * "^'"S that a leak haa bean 

ir: :erar:e'„Tarars~— ^^^^^^ 

mean, tor collecting a plurality ot measuret^ent ft™ ^ * °°"'P««' '"""les a processing 

Plurality ol measutament S^^T^^^Z^^TJ^, 7.'"","""" ' ^""9 '"^ 

.hacontpreased n.tnx tonna, to ..^Z^^IX^S^'^T^X '"'"'^ "'^""^ ^' 

KLaX^rir'n^rgru:^:^^^^ '"I't-^ — • =y«- ™v 

ground storage tank. naergrouna storage tank, an above-ground storage tank, or a partially above- 

tank probe. The dIspensVg'apZa urnJiln "ude a t^^^^^ '^"'^ "^^^ ''^ ^ magnetostrlctlve 

apo™ 

aS:u;:"r=^^ 

disposed at different locations within tSe Tystem each^ ZT h? ' ^"'^"^ ^^'"P^^at"^^ ^^nsors 
P.ur..ofte__ents^^^^^^^^^^^^^^^ 

mS. Fu'r^ZTr^ct^rnTJ:;^^^^^^^^ ! --'^ transpose of the data 

matrix with the transpose of each of rpaSons ^ °' °' " P"^''""ty °^ P^^'tions of the data 

iS. The^:;^^^raTuraT5^^^^^^^^^ ^Ti:. ?:t °' ^^'-'^^ ^'^^-^'-^ 

IVIeasurement data is simultaneousircTlfected fTom tS nT f ^ associated with the system, 

computer to detemiine rchanqe in ih^vo^^^^ °^ measurement apparatus in a fom, readable by a 

data from the plur^ftrof meaTelnt ap^^^^^^^^^^^ T '° " ^'^^^''^y "measurement 

fomiat in a computed meToT and th^^cSeLL m^^^^^ measurement data is stored in a compressed matrix 
monitoring Information. compressed matrix fomiat .s stafsticaiiy analyzed to detemilne operational 

[0037] implementation of the invention may also include the following feature. The method may include estimating 
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an initial value of the volume during the analyzing step. 

[0038] In general, in another aspect, the invention features a method of monitoring a fluid storage and dispensing 
system. The system has measurement apparatus for measuring a volume associated with the system and a plurality 
of temperature sensing devices located at different heights in the system. The volume has a height in the system. A 

5 plurality of volume measurement data is collected from the measurement apparatus in a form readable by a computer. 
The volume measurement data is adjusted based on temperature measurements taken from those of the plurality of 
temperature sensing devices at a height below the height of the volume in the system. The plurality of volume meas- 
urement data is stored In a compressed matrix fomiat in a computer memory, and the compressed matrix fomiat Is 
statistically analyzed to determine operational monitoring information. 

10 [0039] In general, in another aspect, the Invention features a method of determining a volume associated with afluld 
storage and dispensing system. The volume has a height in the system, and the system has measurement apparatus 
for measuring the height. A plurality of height measurement data is collected from the measurement apparatus in a 
form readable by a computer. The plurality of height measurement data is stored In a compressed matrix fomnat in a 
computer memory. Regression analysis is performed using the compressed matrix fonnat to calculate the volume 

IS associated with the system. 

[0040] Implementation of the invention may also include one or more of the following features. The collecting step 
may be performed each time a portion of the volume is dispensed from the system. The collecting step may not be 
performed when the fluid Is being added to the system. 

[0041] In general, in another aspect, the invention features an apparatus for detemnining a volume associated with 
20 a fluid storage and dispensing system. The volume has a.height in the system, and measurement apparatus measures 
the height of the volume. A computer has a processing means for collecting a plurality of height measurement data 
from the measurement apparatus and a memory for storing the plurality of height measurement data in a compressed 
matrix format. The processing means performs regression analysis of the compressed matrix format to determine the 
volume associated with the system. 
25 , [0042] In general, in another aspect, the invention features a method of determining a plurality of volumes, each of 
Vi. the volumes associated with one of a plurality of fluid storage and dispensing systems. Each of the volumes has a 
: height In its associated system, and each of the systems has measurement apparatus for measuring the height for 
each of the volumes. A plurality of height measurement data is collected from the measurement apparatus of each of 
. 4... the plurality of systems in a fonn readable by a computer. The plurality of height measurement data is stored in a 
30 ^ compressed matrix format in a computer memory. Regression analysis Is performed using the compressed matrix 
format to calculate the volumes associated with the systems. 
,7;iv [0043] Other features and advantages of the invention will become apparent from the following detailed description, 
. . and from the claims. 

35 Brief Description of the Drawings 

[0044] 

Fig, 1 is a schematic diagram of a facility Including an underground tank storage system. 
Figs. 2, 3 and 4 are a portion of the Mathcad computer code used to perform the data compression algorithm. 
Figs. 5, 6 and 7 are a block diagram of the steps perfomned during routine operation of the algorithm of the present 
Invention. 

Fig. 8 is a block diagram of the steps performed during the data deletion operation of the algorithm of the present 
invention. 

Fig. 9 is a block diagram of the steps performed during the delivery calculation operation of the algorithm of the 
present invention. 

Fig. 10 is a schematic diagram of a data acquisition and transmission network that may be used in conjunction 
with the present invention. 

Fig. 11 is a schematic diagram of a facility including an above-ground tank storage system. 
Fig. 12 is a schematic diagram of a facility including a partially above-ground tank storage system. 

Description of the Preferred Embodiments 

[0045] The method and apparatus described herein applies to UST's, ASTs or any type of storage tank. The product 
55 stored in the tank may be any fluid, including dry particles that flow in the manner of a fluid, 

[0046] Fig. 1 shows a UST facility 10, illustrated as an automobile service station. Facility 10 includes a series of 
UST's 1 2, 1 4, 1 6 which may store the same or different types of liquid fuel product 1 8. Volumetric tank gauges 20, 22, 
24 in each tank measure the height of product 1 8 in the tank. Submersible pumps 26, 28, 30 in each tank pump product 
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18 to one of dispensing pumps 32, 34 througii piping lines 36, 38,40. Alternately, facility 10 may be an AST facility 
with above-ground tank 1000, as shown in Fig. 11 , or a facility with a partially above-ground tank 1010 as shown in 
Fig. 12. 

[0047] Tank gauges 20, 22. 24 are mounted in tanks 12, 14. 1 6. The tank gauges may consist of or be based on 
magnetostrictive tank probes or other sensing technologies. In the case of magnetostrictive technology, two floats 42 
44 surround each probe, e.g., gauge 20 in tank 12. One float 42 floats on the upper surface of product 18 in tank 12' 
and the other float 44 floats on the interface of product 18 with any water or other foreign material collected at the 
bottom of tank 12. Tank gauge 20 determines the distance between floats 42, 44 to obtain the height of product 18 in 
tank 1 2. Tank gauge 20 also contains temperature sensors 46, 48. 50 spaced along its length to monitorthe temperature 
of prpduct 1 8 at various depth levels. 

[0048] Each of the dispensing pumps 32, 34 consists of a totalizer or flow meter 52. 54 disposed in a housing 56 
58 to measure the volume of product 18 dispensed through hoses 60, 62 and nozzles 64, 66. To operate dispensing 
pump 32, nozzle 64 is removed from housing 56, which actuates dispensing pump 32 and causes product 18 to flow 
through hose 60 due to the pumping action of submersible pumps 26, 28, 30. A value stored in totalizer 52 is incre- 
mented as fuel is dispensed through hose 60. Upon completion of the transaction, nozzle 64 is replaced in housing 
56, thereby turning off dispensing pump 32 and discontinuing the action of submersible pumps 26, 28, 30 and totalizer 
52. 

[0049] Transactions are recorded electronically by software in a sales recording device 71 connected to totalizers 
52, 54 of dispensing pumps 32, 34. Totalizers 52, 54 in dispensing pumps 32, 34 are connected to sales recording 
device 71 by means of communications and power supply wires 78, 80. 

[0050] Sales recording device 71 contains software capable of emulating the functions of a point of sale (POS) 
temninal associated with fuel sales made at facility 10. POS emulation software in sales recording device 71 functions 
on the basis of read only commands to eliminate the possibility of conflict with control commands from a POS terminal 
employed by facility 10. Alternative data acquisition systems can result in destruction of credit card sales records, 
inadvertently shutting down the entire system, and/or causing electrical Interference in the pump links. 
[0051] Tank gauges 20, 22. 24 are connected to a tank monitor 82 by means of communications and power supply 
wires 84. 86, 88 or communicate data through radio frequency transmission. Tank monitor 82 converts raw data ob- 
tained from tank gauges 20, 22, 24 into a fonri usable by a computer 

[0052] A computer 70 contains a processor 72 capable of running various computer software applications and a 
memory 74. Tank monitor 82 and sales recording device 71 are electrically connected to computer 70 to relay totalizer 
values, product height and temperature data to computer 70. Software executable by processor 72 of computer 70 is 
capable of querying tank monitor 82 and sales recording device 71 to obtain measurement data at selected time in- 
tervals. The data is continuously evaluated as it is collected and Is stored in memory 74 of computer 70 for later retrieval, 
and detailed analysis. Alternatively, computer 70 may communicate with a host processor 90 at a remote location The 
continuous evaluations or detailed analysis may then be perfomned by host processor 90, which may be faster or more 
efficient than computer 70. 

[0053] As an example, computer 70 may be a personal computer or any other proprietary microprocessor-based 
unit. Computer 70 may capture data automatically through direct-connect serial interfaces with tank monitor 82 and 
sales recording device 71 , or by manual operator keypad entry. Computer 70 communicates with equipment at facility 
10 through four programmable serial communication ports, such as RS-232 communication ports. 
[0054] Computer 70 may, e.g.. store tank dimensions and product characteristics, and concun-ent time and date data 
along with the measurement data. Computer 70 may be used to produce error and analysis reports as calculated by 
the software. It may also have alarm event-initiated capabilities, such as when a leak is detected in any of the tanks. 
Such a computer system can accommodate facility and customer specific requirements while-maintaining complete 
compatibility with other system components. 

[0055] The SIR method involves reconciling volume data obtained from tank monitor 82 and volume data obtained 
from sales records. Sales transactions may be detected in a number of ways, including an electronic signal emitted 
from totalizers 52, 64, by voltage sensing of control relays on pump dispensers 32, 34, or by observation of product 
removal using tank gauges 20, 22, 24. 

[0056] It is essential that the measurements used to obtain these two types of data are made simultaneously. The 
SIR method of the present Invention collects and analyzes observations of sales volumes and tank volumes which are 
derived simultaneously. Failure to collect both types of data simultaneously would bias estimates derived from separate 
volume measurements. 

[0057] The SIR method properly accounts for the effects of temperature, pressure and specific gravity. In addition, 
product from two or more tanks may be blended, such as to achieve varying petroleum octane levels at pump dispensers 
32, 34. When different fluid products are blended, the tanks are. treated as one unit, and an additional parameter is 
introduced to determine the actual blend percentages. 

[0058] Data concerning the physical characteristics of the tank configurations and the accuracy of the various gauges 
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and metering devices is collected during installation and a set-up phase of operation of facility 10 to create a basis for 
subsequent statistical analysis. Information is then continuously collected so that the statistical analysis of SIR can be 
performed by computer 70 or host processor 90. 

[0059] Several procedures are used either singly or in combination to obtain the volume observations. First, where 
5 the system configuration provides for determining whether hoses and dispensers associated with a given tank are 
active, the system is queried on a minute-by-mlnute basis, or on the basis of another predetermined time interval, to 
determine the status of the dispensers. When all of the dispensers are idle, the values from totalizers 52, 54, the tank 
volumes (i.e. product heights in the tanks) and temperatures are recorded. 

[0060] Second, submersible pumps 26, 28, 30 are checked to detennine on/off status. When It is detennined that 
10 the pumps are turned off, the values from totalizers 52, 54 are read, and tank volumes and temperatures are recorded. 
[0061] Third, software algorithms used by computer 70 detect and measure leads and/or lags between the recording 
of sales events and corresponding gauge and meter readings. When leads or lags are encountered and constitute a 
physical characteristic of the data measurement and recording system, constrained optimization, rather than uncon- 
strained optimization, may be used to determine parameter estimates. Lagrange multipliers are one example of such 
15 a constrained optimization method. 

[0062] The method of the present invention is capable of providing dynamic monitoring of system perfonnance. For 
example, the leak detection function is carried out continuously while normal operations, e.g., removals and deliveries, 
are taking place. To detect leaks dynamically, the software is programmed to detect when sales or delivery events 
occur and to calculate the volumes of product removed or added as a result of such activities. Thus, dynamic testing 
20 does not require that the system be domnant and addresses the entire system from the point of filling to the point of 
dispensing. 

[0063] The SI R method of the present invention also distinguishes between one-time removals and continuous losses 
consistent with leakage. The integrity or leak-free status of the system is not assumed a priori. Instead, the individual 
and unique characteristic pattern induced by each form of error when viewed along the separate dimensions of time, 
25 produQt height and sales volume are used to identify and quantify the errors. The method may also be used to detect 
and quantify undocumented removals, e.g., theft or additions of product. 

[0064]>^^ Further, the overall system is self diagnosing in that it determines from the data the maximum degrees of 
reliability and precision of which a particular operating configuration is capable at any given time, as well as the degree 
of calibration accuracy. 

30 [0065]iy;. In particular, product height in the tanks and temperature are measured continuously at, e.g., one-minute 
interva,ls. Heightandgross volumes are converted to net volumes at, e.g., 60**For1 S^'C, using the algorithms described 
below;>?Sales recorded by the totalizers 52, 54 are extracted and stored in memory 74 at times coincident with readings 
from t^nk gauges 20, 22, 24. If the dispensing system is capable of transmitting a signal indicating whether or not any 
or all individual hoses are active, that infonnation is also stored in memory 74 coincident with taking gauge and meter 

35 ■ readings. 

[0066] The method of the present invention is designed to achieve the maximum accuracy possible within the limi- 
tations imposed by the inherent random and irreducible noise in the various measuring devices incorporated. It utilizes 
multiple measurements over extended time periods to identify and quantify systematic and repeatable effects in the 
instrumentation and thereby correct for such effects using the known physical characteristics of the devices. The system 

40 makes no a priori assumptions as to the accuracy . of the devices used to measure product volume in the tank, to 
measure volumes removed, or as to the accuracy of volumes reported to have been delivered into the system. 
[0067] The resulting volumetric calculations are independent of the physical characteristics of the tank configuration 
and the various measuring devices which may be incorporated in the system. The results do not rely on input entered 
externally by the operator or from diagnostics internal to the measuring devices used. Instead, the output produced by 

45 the software which analyzes the measured data depends only the patterns induced in inventory data produced by the 
tank gauges and measuring devices and, in particular, the cumulative variances that result when the various input 
values are combined. 

[0068] Various error patterns which the measuring devices can induce and the effects of temperature, tank geometry, 
and orientation on cumulative variances are derived from empirical analysis of real-worid inventory data. The system's 
50 software synthesizes the output measurements of the various devices based on known characteristics derived from 
the empirical data. Thus, the software is capable of identifying measurement errors caused by the measuring devices 
and simultaneously compensating for the effects of those errors. 

[0069] Gauges can be systematically inaccurate in two ways. The height of the product in the tank can be incon-ect, 
. and the height to volume conversion algorithms may not reflect accurately the true dimensions of the tank or its orien- 
55 tation in the ground. The latter may be the result of incorrect measurements or an inappropriate conversion algorithm. 
[0070] The presence of such systematic effects and their nature may be established by examining the pattern of 
inventory variances as a function of product height. Errors of this kind induce patterns which repeat themselves as the 
tank is filled and emptied. If the tank length is incorrect, a linear pattern is induced. If product height is in error, a 
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achievable m any one individual observation. Since the leak rate is computed from a series of succS^^e obse^trns 
however, the minimum detectable leak rate can be reduced to any desired magnitude by increasinrthe number o^ 
b~m'eS^^^^^ recorded. Thus, the system can se.e as a final verification for leakage indicalions oSed 

l^im.J^!!'^ conclusion of an initial set up period of data collection including one or more delivery and sales cycles 
dimfn data is analyzed by regression analysis. The initial set-up regression isled trderiv^ tink 

d rnens^oris arid onentation. individual meter calibrations and secular trends. A confidence level value pVcompSed 
at the .01 level of significance to detemiine the minimum leak rate detectable by the system and the ^eskJuatv^riance 
IS computed to provide the current noise level of the system. variance 
[0075] The regression is peri'omied according to the following equation: 



where: 



st,(R,L,T) Volume in gallons derived from the ith gauge reading In inches In a cylindrical tank with or without 

hemispherical end caps with radius R. length L, and tilt over its length of T inches 
^ - 'n't'al inventory In gallons, which is to be estimated. 

= Sales volume recorded on the kth totalizer 

^ riSmated^ ^^'^^ ^''^'''^^ °" ^''^'^^"^ '^""^^^^ ^^'^^ t)e 

= Volume of the yth delivery, 

p^i = Elapsed time since initiation until the ith gauge reading Is recorded. 

- Constant gain or loss in product per unit of time 

rr„r.^:rr I^IS '^^"^^^^ °^ ^-'^^ ^^--^e time mterva. prior to or 



during observation period j. 



■^iJ (l if J 2. i 



[0076] All of the parameters are estimated simultaneously using least square estimation procedures The R and T 
parameters are derived numerically, but the other parameters are derived analytically P™'^^'^"^^^' R and T 
mul^n °^ P«'-«'^etsrs. including the initial inventory, are estimated simultaneously. The initial volume 

must be estimated from all succeeding data, even if the tank is initially empty, otherwise the Initial qauae readrno a^d 
Its conversion to gallons is assigned a credibility not assumed for ail succeeding readings Al^t a X r^SStJ of 
applications, the initial inventory in an already existing and operating system is not accuraterknown^ 
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[0078] Initial inventory estimation Is vital in detemiining the geometry of the tanl<. When tank geometry, tan(< orien- 
tation, ortank product height measurement depart from the values obtained from nominal sources, all gauge and meter 
measurements are affected. It is practically impossible to detect the errors induced in the gauge measurements and 
correct for them unless the estimation of the initial inventory is made coincident with the estimation of the values of the 
other parameters. 

[0079] The estimate of the parameters are based on the totality of the data collected. This means, e.g., that the 
estimate of leak rate Ls is determined from a linear trend Including all of the data collected, not merely at one end of 
the reconciliation period. Likewise, estimates of tank dimensions and orientation are derived from their overall contri- 
bution to reduction in residual variance, as opposed to a sale by sale analysis of tank segments. 
[0080] The volume st((R,L,T) is derived from the product height measurement by multiplying the constant area of 
tank segments of height h {in inches) by tank length L. The volume in gallons of product in a horizontal cylindrical tank 
of radius R is given by: 



[0081] In the case of a tilted tank, the area of the segments vary with position along the length of the tilted tank, and 
the volume is detennined by integrating over the length L. Such integration does not result in a closed fomn because 
the cross sections are not circular, and a numerical integration would severely limit the frequency of observations. 
Instead, in this application the tank is treated as lying horizontally and the product is considered tilted, to derive an 
equivalent volume. This integration yields the closed fpmri: 



The integrand is evaluated between the normalized product heights in inches, hu/R and hi/R, at the lower and higher 
ends of the tilted tank, respectively. It is standard Industry practice to install tanks on an incline to divert water and 
sludge away from the submersible pumps. 

[0082] Tank tilt is identified from the pattern it induces in the record of cumulative variances as a function of product 
height. It is compensated for by fitting the con-ect mathematical form for height to volume conversions in a tilted tank 
to the cumulative variance calculated by the method of least squares. This is done simultaneously with estimation of 
the initial inventory. 

[0083] Tank length L and radius R are established by equating the first partial derivatives of the sum of squared 
cumulative variance with respect to length and radius and detemnining the values which minimize the sum of squared 
variances. Simultaneous estimation of initial inventory is also required when estimating tank length L and radius R. 
[0084] Errors in measurement of the product height h in the tank are characterized by curvilinear patterns induced 
by height to volume conversions in the cumulative variance for a cylindrical container when heights are transposed 
upward or downward. Such errors also are compensated for by minimizing the sum of squared .cumulative variances 
with respect to increments or decrements to measured product height. This estimation also requires simultaneous 
estimation of the initial inventory of the tank. 

[0085] In general, the accuracy of the estimates of the tank dimensions, tank orientation and height measurements 
is confirmed by observing that the cumulative variances of each derived value as a function of nominal product height 
are random and display no systematic influence or effects, 

[0086] Dispenser totalizer calibration is continuously monitored and evaluated by minimizing the sum of squared 
cumulative variances with respect to multiplicative constants associated with individual reported cumulative sales vol- 
umes from all pump dispensers associated with a particular tank system. This eliminates the need for manual verifi- 
cation of meter calibration. 

[0087] In particular, gauge performance Is continuously monitored to Identify gauge malfunctions or degradation in 
gauge performance. Monitoring of gauge performance is independent of diagnostics which are internal to the measuring 
device. Diagnoses of problems are based only on their Impact on the cumulative Inventory variances which are con- 
tinuously monitored by the software. 

[0088] If the gauge fails to record changes in product height when the dispensers register sales, an increase in 
cumulative variances approximately equal to sales volume is observed; this effect can be identified by the monitoring 
software and a warning of gauge malfunction generated to the operator. 




1 3 hu 



^ [ (z - 1)sin-^ (2z - 2^) ^ - 1 (2z - " 5^ 



EP 0 904 526 61 



^!rid rrj' T °^ ''f --^g'^^^ri^S P^'duct height change, but with a time lag after sales are 

nr.^t u!?:.T^ be a feature of normal gauge perfomnance. Such nomial gauge performance is IdentLd by repeated 
to on^ir h H wir""""'"!'^" ^= "^'"^ '=°'^P'«^«^ subsequent return of the cumulative 3ance 

to nomial bounds. When such gauge function is detemiined to be the nomial operating characteristic o rrTartteu^^^^^ 

fr::ceT:r.:rrrZ^^^^^ ^^^^-^^ . introduced mto the software^the^^ise, thV^^uXr 

[0090] Finally, temperatures In the tank are monitored to detect changes that are excessive for the time intervals 

[0091 The software computes actual, rather than nominal, delivered quantities and requires no Input by the svstem 
operator. The operator may choose to input into the system the nominal delivery quantity tncI^K the^^^^^^^ 
wm°th';; ' T ^^'"P^^^^"^^ '^^-"--"t °^ -'^P^nsion of the product aTt'he point ^f Sp The s'^^^^ 
w II then compute overages or shortages between the nominal and actual quantities delivered, as well as l^e oteTes 
or shortages caused by temperature-induced variations in the transport of the product to the acility and n the Subst 
quent mixing of the delivered product with that resident In the tenk ^ 
[0092] Deliveiy Is identified by the software when a positive cumulative variance Is observed which exceeds the 
system noise level and is not succeeded by a return to nomial variance bounds. Delivered quantWes are co^^^^^^^ 
by estimating the volume increases they induce in multiple, successive obsen^ations. The requrd number of suedes 

toteraTceThrr 1' 1T""' '"'''"^"^ ^'"''''^ ^ -"^^^ ^^^^ - wthin a prrdeterned 

nnf^T? t J 7^'^"" °^ '^f P'^^^"' ^^P^'^'^ °' accounting for sales conducted during delivery and for 

noise introduced by post delivery turbulence in the tank. "euvery ana ror 

Efi^' remo^,aiB or additions to the tenk are computed in the same manner Deliveries 

exL?H Jon n °' '"P"^' '^'^''^'^ °f gravity deXer^S^ 

s^ri^XTeir^r:^^^^^^^^^^^^^^ 

f!!.?!'*^ ir^^"^^® ^^^^^'^ "^^"""^'^ ^ continuous linear negative trend in the date which exceeds the 

computed minimum detectable leak rate after all of the various error phenomena described above havTbern identif ed 
and corripensated for. This calculation deals wfth the totality of the date obteined by constan«rmonSfknown 

" T'""'"' ? observations made only when the system is dormant. It is also indepen^^^^^^^^ 
single data reconciliation calculation in that trends throughout all of the date are evaluated eP«n°eni any 

[0095] All calculations concerning volumes are made on the basis of net volumes, according to the following defini- 

Net Volume 

in Tank = Gauge Volume (1 - (t-60)CE) 
where: 

t = Measured temperature in degrees Fahrenheit (if centigrade, the term in parentheses becomes (t-1 5)). 

CE = Coefficient of expansion. 

and 

Net Sales 

Volume = Metered Sale 

(1-(!l^.60)CE) 

where ti and tg are temperatures measured by the tank gauge at the beginning and ending of a sale transaction 

[ZSSnr^ 'fh T'"''' " ''"^ '° ''-"t'^i- requi^d based r'Smai 

infonnation input by the system operator, as follows: uncAicmcji 

GT = Gross gallons on invoice at the originating temninal. 

NT= Netgallons on invoice at the tennlnal. 

XT = Temperature at the terminal. 

CE = Coefficient of expansion. 
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[0096] The program also records: 

tA= Ambient temperature in the tank prior to delivery, 
tF = Temperature in the tanl< at the conclusion of delivery. 

[0097] The following value Is computed: 

ts = Temperature of the product in the delivery vehicle at the facility at the beginning of delivery, 
where: 

NVD = Actual net volume delivered, previously computed. 
NVA = Net volume in the storage tank at the start of delivery. 
NS = Net overage(+) (underage (-)) in delivery. 

= NT-NVD 

GVD = Gross volume delivered. 

= NVD (1 +(tF -60) CE) 
GVS Gross volume in the transport vehicle at the facility prior to delivery. 

= NVD (1 + (ts - 60) CE) 

GSM =. Shrinkage due to mixing in the tank. 

= GVS -GVD 

GVT = Actual gross volume in the transport vehicle at the facility. 

= NVD (1 + (tT - 60) CE) 

GST = Shrinkage during transit to the facility. 

= GVT -GVS 

GOS = Gross overage(+) (underage(-)) adjusted for temperature effects 

= GT - GVD + GST + GSM 

[0098] Calculations of volumes actually delivered are based on multiple observations of the balance of measured 
tank volumes and cumulative sales. This method requires frequent simultaneous observations of sales and in-tank 
volumes (i.e. product heights) and temperatures, 

[0099] The volume of product in a tank is derived by measuring the height of the product and. using the geometry of 
the tank which is assumed to be known, to compute the corresponding volume. In many instances, tank dimensions 
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vary substantially from assumed design dimensions. Regulatory specifications permit up to 10% variation in ienoth 
and diameter of cylindrical tanks. ^ 
[01 00] Tank orientation can also cause complications in the calculations. The volume corresponding to a measured 
height varies substantially when the tank is tilted away from horizontal or rolled away from vertical. 
[01 01 ] Further, tanks may also fail to confonn to a known geometry either through faulty manufacture or installation 
or may suffer significant defonnation during the course of operations. For example, many fiberglass tanks saq or bend 
along their length. 

[0102] In addition, installed tanks are typically inaccessible, and difficultto measure. Thus, it is necessary to confirm 
the accuracy of height to volume conversions from generated inventory data and to Identify and correct discrepancies 
where they exist. 

[01 03] The foregoing problems are compounded when two or more tanks are manifolded together. Manifolded tanks 
are joined together by piping systems and serve common dispensers. Thus, sales quantities from manifolded tanks 
constitute withdrawals from all tanks In the manifolded system, but not necessarily in equal quantities. Product heights 
typically vary from tank to tank, but tank geometries, dimensions and orientation may also vary so that a procedure for 
correcting height to volume conversion errors for a single tank will not apply 

[01 04] The different factors which Influence inventory data manifest themselves In distinct ways which facilitate their 
Identification and correction. These factors are most easily identified by examination of their effects on cumulative 
departures of actual measured inventory from a theoretical or book value when viewed across a variety of dimensions 
In particular, one-time undocumented physical additions or removals of product, e.g. over or under deliveries and 
pilferage, are evidenced by an addition or subtraction of a constant quantity from the cumulative variance at the time 
of occun-ence and all subsequent observations. Continuous loss of product accumulating overtime, e.g. leakage is 
evidenced by a loss trend over time. Continuous loss of product varying proportionally with sales value, such a line 
leak or meter miscalibration. may be determined by identifying a constant negative trend that is cumulative only over 
periods where delivery lines are pressurized. 

[0105] A pattern of gains or losses, or both, recurring cyclically as the tank is successively filled and emptied with 
no long term gam or loss of product, Is the pattem associated with height to volume conversion error. The pattern is 
cyclical because the error source Is identical in each cycle as the tank Is filled or emptied. It is distinguishable from the 
other patterns in that it retraces the same path without the translation which would occur If physical loss or aain of 
product were taking place. aam 
[0106] This problem is most readily diagnosed by analyzing cumulative variance as a function of product height If 
the vanances are random with no evidence of systematic effects, height to volume conversions may be assumed to 
be correct. If not, the fomi of the induced pattem indicates the nature of the conversion error. Thus, an error in tank 
length induces a linear pattern, an enror in tank tilt induces a sinusoidal pattem, and a constant error In tank height 
measurement Induces an arclike pattern. When other sources of loss or gain are present, the conversion error patterns 
remain, but are translated In each succeeding filling/emptying cycle to reflect the physical loss of product which has 
occun-ed dunng that cycle. Thus, confusion between conversion en-ors and other effects can be eliminated 
[0107] Sales readings and product height measurements must be made simultaneously. Since the number of ob- 
servations in any one sales cycle is typically too few to generate a conversion table of sufficient detail to be of practical 
use. subsequent sales cycles and their corresponding deliveries must be Incorporated. If. however deliveries are 
unmetered and are used to approximate the volume (as is the standard industry practice), significant inconsistencies 
are introduced. If an overage or shortage occurs during delivery, then all subsequent sales volumes correspond to tank 
cross sections which have been shifted upward or downward from their predecessors. Averaging or statistical treatment 
cannot overcome this deficiency since there is no means of knowing without metering whether, by how much and in 
what direction the data has been shifted. . 
[0108] The procedure of the present invention may include detemiining if height to volume conversion error is a 
problem. If the error is a problem, then the system must determine the nature of the problem, e.g. tank dimensions 
tank onentation, height measurement or unknown tank geometry, and whether the conversion problem is compounded 
by other gams and losses. If leakage is suspected, an on-site leak detection investigation is undertaken. In no leakage 
IS indicated, and one or all of tank dimensions, tank orientation and height measurement are problems, new conversion 
factors are calculated and confinned using the diagnostic procedures described herein. 

[0109] If unknown tank geometry or manifolded systems are encountered, the exact current percentage of metered 
sales actually dispensed from each dispenser is determined by physical measurement. A high order polynomial using 
a vanab e of measured product height is used to convert height to volume. The parameters of the polynomial are 
derived from the differences between measured product height con-esponding to the beginning and ending of sales 
events which do not overiap deliveries. 

[01 10] For a single tank, actual dispensed quantities are regressed using a polynomial based on the differences in 
measured product height before and after individual sales, subject to the constraint that when the polynomial is eval- 
uated at a height equal to tank diameter, the result is the total tank volume. Observations which include delivery events 
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ASale, = (h^^ - h,) + 33(11^ - ,) + 



Voi = a.d + a9d^ + ... + a„d" 



[0111] A fifth order polynomial has proven adequate in most cases. Residual analysis may be used to determine 
adequacy of the polynomial in the presence of severe tank distortions, and higher order polynomials may be introduced 
as necessary. The number of observations required Is detemnined by estimating a confidence bound around the re- 
sulting polynomial with a width adequate for the desired resolution. Thus, 



ASaiej = Actual dispensed volume In period /. 

hj = Product height upon conclusion of ASale,. 

20 hj.^ = Product height prior to commencement of ASalei and after completion of ASalej.v 

d= Diameter of tank. 

Vol = Total volume of tank. 

The converted volume for height h is then given by: 

■'"If- 

. Vol (h) = a^h + agh^ + ... + a^h" 

The orriission of a constant term in the regression implies that 

Vol (h) = 0 when h = 0 
This ensures that the polynomial derived from the height differences is well defined. 

[0112] For manifolded systems, actual sales are regressed simultaneously on individual polynomials based on the 
various height differences in the several tanks which correspond to a particular sales volume, subject to the constraint 
that each polynomial evaluated at the corresponding tank diameter yields the total volume of that tank. 

ASalOj = a.,^{h^ii - h^) + a^^(h^. ,.1^ - h^ 

+ "' + a„i(h"Mi-h"M) 

45 2 2 

h 12) 

. + ... + an2(h" ■ 12) ~+ - 

.+ aim(hMm - ^^lm)> a2^(h^ Mm - Im) 
+ + an^(h" Mm-h"|m) 

Where: 

ASale, = Actual Sales volume in period /. 

h|.^j = Height of product in tank / after completion of Asalej.^ and. prior to commencing Asale,. 
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y= 1,2, ... m 

h|j = Height of product in tank y after completion of ASale,. 

m = Number of tanl<s manifolded. 

Volume conversion for the m measured heights, hi, hg, .. . in the total system is: 



Vol (hi, hj, . . . h„) = z r 



where: 

h| = Height of product measured in the Ah tank in the manifold. 

[0113] Deliven^ inaccuracies have no impact on this calculation since all observations made during deliveHes are 
discarded. Height changes are related only to the corresponding volumes dispensed 

[0114] Prior detemnination of actual quantities dispensed, as opposed to metered quantities, ensures that the only 
reniain.ng source of error is random measurement error. Regression is designed to accommodate random error of this 
kind and to facilitate inferences when errors are present. 

[0115] With respectto temperature, the temperature of product delivered into a tank system almost invariably differs 
rom the emperature of the product already In the tank. Its addition has the effect of expanding or contracting the 
volume of the combined product. This change in volume can create the appearance of incorrect dimensions of the 
rn-!«, ,M '^o'^version. appear as leakage where none exists, or it can mask the existence of actual leakage 
[0116] It Is therefore preferable, and frequently essential, that all volumes, sales, deliveries and product in storage 
be converted to a common temperature prior to analysis. Typically SCF (15-C) is chosen as the standard. The con- 
version IS accomplished as follows: 

Net Volume = Gross Volume (1 - {t-60)CE) 
where: 

t = IVIeasured product temperature in degrees Fahrenheit. 

CE = Coefficient of expansion. 

As above, all calculations are in net gallons of product. 

[0117] A complication to the calculation may occur If the tank gauges 20, 22. 24 used to measure product volume 
are designed for static or domiant mode tank testing. Such tank gauges detect leakage when the tank Is taken out of 
service^ n this case, product volume changes due to temperature changes during the couise of a test must be ac- 
countou Tor. 

^"^^T',^ ^ • '^^P^^^t^i-e sensors 46, 48, 50 are located at different heights in tank 1 2. If the 

drnnn?H^frnl"th Z ^ '^"'P^'^'^^ ^ensor, the Corresponding weighted temperature measurement is 

dropped from the average temperature calculation, and a temperature jump and corresponding volume change may 
be observed when the net volume is calculated using the new weighted average of temperatures. If uncorrected sTc!! 
Xr^e cinTerslons '''"'"'^^ °^ "^'^ '^^^ °' generation of height to 

thrf!LwIn^g%^SL°' '"^^^"^ invention may be used to overcome these temperature related problems. Using 

NDBn = The net cumulative variance in the Inventory data at obsen/atlon N. 
then, 
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NDBm « a(l (to - 60) CE) 

N 

Sai(l - (ti - 60) CE) 

i-1 

-V^d-Ct^ - 60)CE) 



where: 



a = Gross initial inventory, 

to = Tennperature of initial product volume. 

15 t' = Tennperature of product at observation /. 

Saj = Gross volume sold in period /. 

V|vj = Measured gross volume in tank at period N. 

CE = Coefficient of expansion. 

20 Absent random error or leakage, and. assuming no deliveries of product, then 
NDBn = 0 and 



30 ; 



a(l - (to • 60) CE) 



N 

r Sa. (1 - (tj^ - 60) CE) 

i-1 



(1 - (tjj - 60) CE) 

Therefore, if a temperature jump to temperature t* occurs at an observation N + 1 , then 

NDB^i^i=VN(1-(tN-60)CE) 
-Sa^^^ (1 -(t*-60)CE) 
■ Vr4+1 - (t* - 60)CE) 
= Vn(1 - (tN - 60)CE) 
- SaN^-i (1 - r - 60)CE) 
-(VN-SaN^,)(1 -(t*-60)CE) 
= VN(tMN)CE 



[0120] When this final quantity NDB^+i is added to the volume where the transition occurs between temperature 
sensors, and all subsequent volumes, the effect of the transition is eliminated, and analysis proceeds as it would where 
individual temperature readings are available. 

[0121] A large number of variables must be estimated by the software to implement the SIR system of the present 
55 invention. For example, as many as forty hoses and independent totalizers per tank system, as well as deliveries 
numbering four or more per day must be accommodated. Thus, a very large volume of data must be accumulated, 
encompassing a substantial spread of sales volumes from each totalizer for both the set-up analysis and subsequent 
routine monitoring. To accommodate this volume of data within current or conceivablefuture practical computer memory 
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suS21'n?l!llL^?°''l^'" implemented by the software utilizes a matrix formulation which invokes the property of a 
sufficient statistic to reduce the memory requirement. H<v^H<^'iy » 

[0122] The calculations used to determine the various error, loss trend and delivery estimates have the fonn: 

B = (x"^x)"''x*^y 

Hcp- (y-xB)^(y-xB) 
m+1 



= (x'^x)*'' MSE 



where: 



B = Column vector of m parameters to be estimated. 



X = 



Matrix of parameter coefficients. 



Column vector of Independent variables. 



MSE = Mean squared error. 

S2 r= Variance covariance matrix of parameter estimates. 

in'Inml contained in vector y comprise tanic gauge readings. The entries in matrix x are measured sales 

vo umes time, and other constant values. The parameters of vector B which are to be evaluated include the Initial 
volume of the system and subsequent volume changes, including delivery amounts 

vr.tnLrL^''^'^^'!'^ observations are recorded every minute, as many as 1440 rows in the x matrix and the y 

oeS Tnr^: in? H H ^° "'^^"^"'"^^ ^'^'^ ^^^^ ---^ an extended 

period of time. Instead, data compression techniques are applied so that only a manageable amount of data need be 

?i<f f 1 • The algorithm utilizes the property that if an n x m matrix A is partitioned into two submatrlces, B and C. where 
B IS an I X m matrix and C is a j x m matrix, such that i + j =: n, then wnt^re 



c"^c = a"^A -I- b"^b 



[0126] For example, 



if C = 




then 




and B 



2 1 
1 2 
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[0127] At the conclusion of each 24 hour or other period, only x^x and x^y are computed and stored. The nnatrix x 
has the form of a square n x n matrix. Further, the aggregates of observations for different periods are additive, since 
two square matrices having n x n dimensions may be added. Thus the total data storage requirement for each period 
id determined only by the square of the number of parameters of interest. 

[0128] The system Is able to accomodate virtually unlimited numbers of observations by this method of data com- 
pression. Without this capability, the system would not have the storage capacity to accurately and simultaneously 
estimate the numbers of parameters which are required to perform a statistically significant calculation. This data 
compression method also allows for processing the data at the facility or for transmitting the data to a host computer 
for periodic analysis. Figs. 2, 3 and 4 show the Mathcad computer code used to perfomn the data compression algorithm. 
[0129],: Furthermore, (x'^x)-"'x"''y is a complete and sufficient statistic for B, No statistically useful information is lost in 
the compression. The overall procedure is, therefore, unlimited by memory. The only limitation remaining is the precision 
available in the computer system used. 

'[0130]. . The software perfonns SIR analysis, including inventory estimation and leak detection, using the above equa- 
tion in the following form: 

X= [Iranqe. (Trenge). ("D (S range 



where: 

range = 1 ... (number of observations) 
meters = 1 ... (number of dispensers) 

and 

Column of Vs. 
Cumulative time in minutes. 

Cumulative sales for an individual dispenser in gallons. 
Cumulative deliveries. 
Tank stick reading in gallons. 

[0131] To estimate the initial inventory, the matrix x includes a column of unitary values. To estimate loss trends; the 
matrix x includes a column containing cumulative times of measurement and cumulative sales. The values of B, MSE 
and S2 are then calculated, producing the following result for the vector B: 



^ range ~ 

''"range ~ 

/c 's<meters> 
V'^^range^ 

^^rang© = 

Stk^ange ~ 
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^1= Estimated initial Inventory. 

^2 = Loss trend. 

S2+meters = Individual meter error. 



S2 is the variance covariance matrix of the parameter estimates THmq <5 - /c2 m/2 • *u * 
[01 33] The software perf omis delivery calculations using the equation in the following fomi: 



^ = Hrange. (Trange) (-1), (Srange) (-1). (Drange)] 

y - (Sti<,.^p,gg) 
where: 

range = 1 ... (number of records) 
and 

"'range = Column of 1*s. 

Cumulative time in minutes, 
^range Cumulative sales in gallons. 

[0134] Tho .alues of e, MSE ana are th.n caloulaed, producing ,he following result to, the vector 

= Estimated initial Inventory. 
Bg = Loss trend. 
B3 = Meter error. 
64 = Estimated delivery amount. 

[0135] 32 is the variance covariance matrix of the estimates Thus S - ^<52 m/2 io th« ^ ^ -1 . . 

^24 %8'"1?2T3^and^T^T °" '^''t °' ""'^"'"''"^ ^'^^ °f CounteM in steps 114 120 

124. 128, 132, 136 and 140. Depending on the cumulative variance and tiie value of Co..ntf.ri th« \ 

upon updating the value of Counterl and other computational variables (I e index indivlnH ci^lf i ^ 
collected data is deleted (steps 122 and 134). "^n^^'SMi e. index, sign index andsign). In some cases, 
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[0140] Upon computing the rate loss at step 144, the program reads new data at step 146 if the rate loss is not 
greater than or equal to, e.g., .2 gallon per hour, otherwise It computes the trend of the data at step 148. If at step 150 
it is detemnined that the trend is greater than .2 gallon per hour, a warning is issued at step 156. In either case, the 
software continues to read and analyze the data at steps 152, 154, 158 and 160 until the last observation. 
[0141] The operation of deleting data 1 70 is shown in detail in Fig. 8. After performing similar analyses at steps 1 72, 
1 74, 1 76, 1 80 and 1 84, using the indices and the values of the calculated standard deviations as in the routine operation 
procedure described above, the values of the counters are updated and new data is read at steps 178, 1 82, 1 86 and 
1 88. Data is deleted in accordance with steps 1 78, 1 86 and 188. 

[0142] Finally, Fig. 9 shows the delivery calculation 190 in detail. After determining that the cumulative variance is 
greater than a predetermined value (three standard deviations) at step 1 92, the program detemiines whether the var- 
iance is greater than, e.g., 100 gallons per minute (step 194). If so, the delivery is recorded and the amount of delivered 
is determined at steps 202 through 220. 

[0143] if there is a delivery in progress (step 202), data is read until a negative observation to observation variance 
is observed (step 204). The variance is monitored until the turiDulence in the tank subsides (step 206). Thirty observa- 
tions are read, (step 208), and all observations from 15 minutes before the delivery until the end of the turbulence 
observations are deleted (step 210). An indicator variable is introduced with the turbulence observation, from which 
regression commences (step 212). The confidence bound on the indicator is computed (step 214), if the confidence 
bound is within a predetermined tolerance, the volume of the delivery is reported within the confidence bounds (step 
220); otherwise, additional observations are added, and the confidence bound is recomputed (step 218). 
[0144] If the variance between data measurements is less than 100 gallons per minute, the software detennines 
whether the gauges are inoperative and reports them as being inoperative (step 198), or proceeds as in the routine 
operation procedure according to step 200 (in which there is a negative variance) depending on whether the observation 
exceeds a predetermined value (within one standard deviation) at step 1 96. 

[0145] in general, if observed variances are within three standard deviations or other predetennined tolerance of the 
expected value, the data Is stored for future analysis, VVhen cumulative variance exceeds three standard deviations 
or other predetennined tolerance, different software programs are executed depending on the nature and magnitude 
of the. departure. 

[0146] If within. ten (or other predetennined) successive observations after the initial departure, the cumulative var- 
iance returns to within the tolerance range, all data from and including the initial departure and prior to the initial ob- 
servation are deleted. The time extent and number of observations involved is recorded and stored for, e.g., a daily 
gauge perfonnance report. 

[0147]> if all ten (or other predetermined) successive observations remain outside the tolerance bound and the cu- 
mulative variances are of the same sign, a new trend line is initiated at the point of initial departure. After ten (or other 
predetermined) additional observations, a third trend line is initiated. If the increment to the overall trend estimated 
from the most recent obsen/ations is not significant, the most recent data is consolidated with the previous data and 
the process is repeated until such time, if ever, that the current trend increment is significant. 

[0148] If the departure is positive, the system checks whether the product is being dispensed and whether the gauge 
height fails to decrease, reflecting removal from tank. If so, the tank gauge is reported to be inoperative. 
[0149] If the gauge height is increasing, monitoring is continued as above until the most recent trend line returns to 
its original slope. Minute to minute variances are monitored to detect turbulence until the gauge values again return to 
within tolerance. All observations which occurred in the fifteen minutes prior to first positive departure until the end of 
post delivery turiDulence are deleted. An indicator variable is introduced at the first observation after post delivery 
turbulence. The system collects thirty additional observations and performs the regression from the beginning of the 
period to determine the volume delivered. The volume delivered is then reported. 

[0150] If the departure is negative, the system proceeds as with delivery. If successive slope increments fail to show 
a return to the original slope, Indicating continuing loss of product for a predetermined period, typically one hour, and 
slope exceeds .2 gallon per hour, the system reports a warning that there is a continuous loss of product. If the loss 
rate is less than .2 gallon per hour but greater than the minimum detectable leak, the system continues to monitor and 
recalculate the parameters, to be included in a daily operational report, if the incremental trend line shows a retum to 
the original trend, the system proceeds as with delivery, introduces an indicator variable, deletes data as necessary, 
and performs the regression to determine the volume of product removed. The system reports a one-time removal of 
product. 

[0151] Referring to Fig. 10, the invention incorporates a data acquisition and transmission networt< (DAT Network) 
300 to completely automate the process of obtaining, capturing, transferring and processing product inventory data 
for use in product management, delivery scheduling and environmental compliance practices. DAT networi< 300 in- 
cludes on-site processors 302, 304 at the facilities 306, 308 where the tanks are located, a customer host processor 
31 0 and a central host processor 31 2. DAT networi< 300 links multiple remote facilities 306, 308 to central host processor 
312, which performs the SIR analysis. The link may be accomplished indirectly through customer host processor 31 0, 
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Which itself is connected to a plurality of remote facilities 306, 308 Each of these orocessor ^.temon^ ^ . 

independently operatingsoftwareandhardwaresystemswhichformSb^ 

Which transmit information electronically via the telephone network 314 using standard dialTvo ^i gr^l Shone 

[.n?emroun''Ln7ir ^'"^f" ^ °^ '"^^^rated management for the widest possible variation of 

underground and above-ground fuel storage, movement and measurement systems. On-site processors 302 304 «rJ 
capable of obtaining Infomiation from any electronic or mechanical control system enabHnq DA^n^ro^ f 00 

Sis S'To?; *° ^-"'ty While presenting the InSSatioLaZ^i at^ 

aciNWes 306 308 to customer host processor 31 0 or central host processor 31 2 in a unifbmi fomiat 
[0153] On-line processors 302. 304 obtain and capture product Inventory data through the use of proorletarv inter 
aces wrth external systems in use at remote facility 306, 308. such as tank gauges and sates reco/dKvS ?^^^ 
neLoT^li!? T"'"' \nionn^tlon dally, weekly or monthly through the pub^ s^Sed IZhone 

scJ!d in « ^7^''°' ° °' P^'^^ssor 312 for use in inventory management d^S 

Sfr"*!. '? P^'^'^"'^'" f ""-^'^^ processors 302, 304 typically is equipped with an alphanumeric keypad a char- 

Ja^ moSmLTr?"'': ^* pons, an interna" auto-dial/autoansw'r (^^^ 

«ni ™. , ? , (^^-^ connection to a printer). The keypad and display allow for operator confruraSon 

and manual entry of sales, deliver, and tank level data. Use of an AD/AA 2400 baud modem alLs Lttil on^sSe 

7:ZZ ' f ^ *° ^" "'''^'"^ ""'^^ ^"'^^ ""^ est^Usi^mg communicarn widows o min^^^^^^ 

attemptedsrmultaneous use. Each oftheprogrammableserialcommunication ports is independent, fullyprir^^^^^^ 
and governed by options selected atthefacility or off-site through modem access. Finally, on-site proceLoJs 302 304 

f^m 1, J . f/ P^'^^^^^^ 31 0, Which is capable of receiving, storing and processing inforr^atlon 

nlZnf^T!r:T ^r"" ^"^'^'^^ ^^^^ management of a remote tank population fro^a S^poim 

of contact. A database of information created by customer host processor 31 0 is the basis for all higher level oroduc 
management functions perfomied by DAT network 300. The database is also the basis for the env ?onmS cor^o^ 
ance analysis performed by central host processor 312. environmental compli- 

fnlhrLt"^^ *^°«» P^°=essor312, Which is capable of receiving, storing and processing the infomnation 

^ the database created by customer host processor 31 0 for product management enables DAT netJoSo to aThteve 
m^rnurn results by^utihzing the databasefor environmental compliance without additional remrfadm? infon^a^^^^ 
an^ivr^r ?f °"- 312 is capable of transmitting a resulting database of the environmS 

analysis back to customer host processor 31 0 for printing and other customer record keeping requirements 
[015^ The processor elements of DAT network 300 may exhibit other useful operational characteristics To orevent 

s™2? ''j''' """"""^ " ''"""^ ^'"'''^ '^'^'-"P data'transfer fun« onTis r^ulJd'u^^^^^^ 

secured access the baud rate, parity, stop bit parameters and communication protocol are determined at Tnv of on 
site processors 302. 304, customer host processor 31 0 or central host processor 312 ''^'^'^'"^^ °f °"- 

[0158] Another function of DAT network 300 Is to monitor tank contents generally. DAT networi< 300 can be 
grammed to activate, e.g.. an audible and visual alami if the water level in The tank is too hir^g greater than 2 

=r IS e^^nS. — °' '-^^ - --^ 

St1!?l y*"^ ^"«'y='^ '"ay also be combined with additional inputs of holding costs, reorder costs transoor- 

tat on costs and penalty costs for running out of stock. The system can include optimal inventory Sor^hms to dSl 

echT °'nt'''T"''"r ^^°^^^^P°'"^^ °P^'^-' ^-"-^ t-'^X routing. Further, the syster,^^ n^oSe mu" 
echelon, optimal inventory procedures to accommodate combined wholesale and retail operations sucHs wTca - 
culus-based optimization and linear, nonlinear and dynamic programming. operations, such as with cal- 

[0161] Other embodiments are within the scope of the claims. 
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Claims 



A method of monitoring a fluid storage and diispensing system, said system comprising measurement apparatus 
for measuring a volume associated with said system, said method comprising; 



collecting a plurality of measurement data from said measurement apparatus in a form readable by a computer; 
storing said plurality of measurement data in a compressed matrix format in a computer memory; and 
statistically analysing said compressed matrix fonnat to determine operational monitoring information. 

2. The method of Claim 1 further comprising 

a) determining the presence of operational defects in said system; or 

b) monitoring the accuracy of said measurement apparatus; or 

c) determining the volume of fluid in said system; or 

d) determining whether a quantity of fluid removed from said system is caused by a leak in the system, and 
preferably, delivering a warning that a leak has been detected in said system; or 

e) querying said measurement apparatus under the control of said computer; or 

f) measuring the temperature of the fluid; and calculating said volume based on the temperature of said volume 
in said system and, preferably, wherein said calculating Includes determining a correction value based on a 
weighted average of the temperature of said volume measured at a plurality of locations within said system; or 

g) transmitting said measurement data to a host processor to perfonn said statistically analysing steps; or 

h) transmitting said compressed matrix fonnat to a host computer to perform said statistically analysing step. 

3. The method of Claim 1 wherein 

■:' . a) said statistically analysing step includes calculating error data resulting from measurement apparatus; or 
V;' b) said collecting step is performed while said system is operating and, preferably, said collecting step is 
j.'t performed continuously at periodic intervals, and more preferably said collecting step is perfonn ed automat- 
'2. icaliy under the control of said computer; or 
. :^v c) said system comprises a plurality of tanks; or 

d) said system comprises an underground storage tank; or 
, e) said system comprises an above-ground storage tank; or 
.'^.^ f) said system comprises a partially above-ground storage tank; or 

g) said storing step comprises generating said compressed matrix format as a product of a data matrix and 
the transpose of said data matrix and, preferably, wherein said product is formed by addition of partial products 
of each of a plurality of partitions of said data matrix with the transpose of each said partition; or 

h) said measurement apparatus includes a volumetric gauge, a dispensing apparatus and a sales recording 
device; or 

t) said collecting step comprises estimating an initial volume of fluid in said system. 

4. A method as claimed in Claim 1 wherein said measurement apparatus comprises a plurality of measurement 
apparatus, said method further comprising the steps of 

a) simultaneously collecting measurement data from said plurality of measurement apparatus in a fonn read- 
able by a computer to determine a change in said volume; 

b) repeating said collecting step to obtain a plurality of said measurement data from said plurality of measure- 
ment apparatus; and preferably 

c) estimating an Initial value of said volume during said analysing step. 

5. A method as claimed in Claim 1 wherein said system further comprises a plurality of temperature sensing devices 
located at different heights in said system, said volume having a height in said system, said method further com- 
prising: 

adjusting said volume measurement data based on temperature measurements taken from those of said plu- 
rality of temperature sensing devices at a height below the height of the volume in said system. 

6. A method as claimed in Claim 1 of determining a volume associated with said fluid storage and dispensing system, 
said volume having a height in said system, said measurement apparatus being for measuring said height, said 



r. 
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statistical analysis being regression analysis using said compressed matrix fomiatto.calculate said volume asso- 
ciated with said system; and preferably said method further characterised In that 

mo're prefer^lj^ ''^^ psrfomied each time a portion of said volume is dispensed from said system; and 
b) said collecting step is not performed when said fluid Is being added to said system. 
A rnethod as claimed In Claim 6 wherein said volume is a plurality of volumes, each of said volumes associated 

Tster: ^eS'^:!?''?"'''"'''^^^^^ 

vofumes compnsing measurement apparatus for measuring said height for each of said 

An apparatus for monitoring a fluid storage and dispensing system, said apparatus comprising: 
measurement apparatus for measuring a volume associated with said system- and 

a computer comprising a processing means for coiiecting a plurality of measurement data from said meas- 
urement apparatus and a memory for storing said plurality of measurement data In a compressed matrixformat; 

.n.J,nnT ^^l^ ^^^f '"^ "^^^"^ P^'^°'^= «t««stical analysis of said compressed matrix forrriat to detemilne 

operational monitoring information. 

■ «n"«LT'^*"^ Claimed in Claim 8 wherein said volume has a height in said system, wherein said statistical 
Sd syieT'''^' ^"^ ^ '^'^ °P"^^''°"^' infomiatlon is said volume associated S 

0. The apparatus of Claim 8 wherein 

a) said system comprises a plurality of tanks; or 

b) said system comprises an underground storage tanl<; or 

c) said system comprises an above-ground storage tank; or 

d) said system comprises a partially above-ground storage tank- or 

e) said measurement data is collected simultaneously from said'varlous measurement apparatus- or 

f) said measurement data obtained from said measurement apparatus are transmitted. to said computer or 
analysis comprises a host processor remote from said computer to perfomi said statistical 

. The apparatus of Claim 8 wherein said measurement apparatus includes a volumetric gauge, a dispensing appa- 
ratus and a sales recording device and preferably . a ui&pensing appa 

a) said dispensing apparatus comprises a totallser; or 

b) said sales recording device simulates operation of a point of sale tenninai- or 

tltrii^ive ta^IIk proS"''' ^ ^^'"^ ^^''^ « '"^a"^- 

. The apparatus of Claim 8 further comprising a temperature sensor disposed within said system in contact with 
said volume for obtaining a temperature measurement of said volume and. preferably, wherein TaS terJo^aru e 
measurement is used to calculate said volume. ' r- nerein saia Temperature 

. The apparatus of Claim 8 wherein 

a) the apparatus further comprises a plurality of temperature sensors disposed at different locations within 
said system, each of said senso,^ being In contact with said volume for obtaining a plurality of tempeSu^ 
measurements of said volume and. preferably, wherein said temperature measurements areled to caSaS 
saiu voiumG; or 

b) said compressed matrix fomiat is generated as a product of a data matrix and the transpose of said data 

n^fonc "f H H t ^' °' P^°ducts Of each of a plurality of 

partitions of said data matrix with the transpose of each said partition. 
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PatentansprUche 

1. Verfahren zum Oberwachen eines Flussigkeitsspeicher- und Ausgabesystems, wobei das System eine Messvor- 
richtung zum Messen eines dem System zugeordneten Volumens aufweist, wobel bei dem Verfahren: 

eine IVlehrzah! von Messdaten aus der Messvorrichtung in etner computerlesbaren Form gesammelt wird; 

die Mehrzahl von IVlessdaten in einem l<omprimierten Matrixformat in einem Computerspeiciier gespeichert 
wird; und 

das l<ompr(mierte Matrixformat statistiscli analysiert wird, um BetriebsubenA^acliungsinfomnationen zu bestim- 
men. 

2. Verfaliren nach Anspruch 1 , bei dem weiter: 

a) das Vorhandensein von Betriebsdefelcten in dem System festgestellt wird, oder 

b) die Genauigiceit der Messvonriclitung ubenrt^acht wird, oder 

c) das Voiumen der Flussigl^eit in dem System bestimmt wird, oder 

d) festgestellt wird, ob eine aus dem System entfemte Menge an Flussigkeit durch ein Leek in dem System 
verursacht ist, und vorzugsweise eine Warnung geliefert wird, dass ein Leek in dem System festgestellt worden 
ist, Oder 

. e) die Messvorrichtung unter Kontrolle des Computers abgefragt wird. oder 

f) die Temperatur der Flussigkeit gemessen wird und das Voiumen auf Basis der Temperatur des Volumens 
: in dem System bereehnet wird und die Berechnung vorzugsweise die Bestimmung eines Korrekturwertes auf 
. Grundlage eines gewichteten Mitteis der Temperatur des gemessenen Volumens an einer Vieizahl von Orten 
innerhaib des Systems beinhaltet, oder 

. ;v,;rv g) die Messdaten zu einem Host-Prozessor ubertragen werden, um die statistischen Analyseschritte auszu- 
- fuhren, oder 

h) das komprimierte Matrixformat zu einem Host-Computer ubertragen wird.. um die statistischen Analyse- 
schritte auszufuhren. 

3. Verfahren nach Anspruch 1, wobei 

a) der statistische Analyseschritt das Berechnen von aus der Messvorrichtung resultierenden. Fehlerdaten 
umfasst, Oder 

b) der Sammeischhtt durchgefuhrt wird, wahrend das System im Betrieb ist, und vorzugsweise der Sammel- 
schritt kontinuierlich in periodischen Abstanden durchgefuhrt wird und besonders bevorzugt der Sammelschritt 
automatisch unter Kontrolle des Computers durchgefuhrt wird, oder 

c) das System eine Mehrzahl von Tanks aufweist, oder 

d) das System einen unterirdischen Speichertankiaufweist, Oder 

e) das System einen oberirdischen Spetchertank aufweist, oder 

f) das System einen teilweise oberirdischen Speicliertank aufweist, oder 

g) der Speicherschritt das Erzeugen des komprimierten Matrixformats als Produkt einer Datenmatrix mit der - 
Transponierten der Datenmatrix umfasst, und .wobei vorzugsweise das Produkt durch Addition von Teilpro-. 
dukten von jeder aus einer Mehrzahl von Teilmatrizen der Datenmatrix mit der Transponierten Jeder Teiimatrix 
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gebildet wird, Oder 

I) der S.mmetechritt Oas Schatzen sines Anteng»,olumens an Flilsslgtelt in dem System umfasst. 

ire"7srzrrrursrr.Xmr^"'°""'^^^^ 

r~;'u" ™ .on Mess.c*nt.n. 

c) Abschatzen eines Anfangswertes des Volumens wahrend des Analyseschritts. 

CSiTnn T ^"'^"^f-^^^^^^t^" Gmndlage der Temperaturmessungen, die von denjenigen aus der 

:l?dlT~r^^^^^^ - Svstem ausgegeben 

b) der Sannmelschritt nicht ausgefuhrt wird. wenn dem System Fiusslgkeit hinzugefugt wird. 

. Verfahren nach Anspruch 6. wobei das Voiunnen durch eine iVIehrzahl von Volumen gebildet wird wobei led«. 
Vo uZ 7nTJT Flussigkeitsspelcher- und Ausgabesystemen zugeSS eT^t' Tobe edes 

. vorrichtung zum Uberwaohen eines Flussigkeitsspeicher- und Ausgabesystems. wobei die Vorrichtung aufweist: 
eine IVI ess vorrichtung zum IVIessen eines dem System zugehorigen Volumens, und 

einen Computer mit einer Prozessoreinrlchtung zum Sammein einer Mehrzahi von Messdaten aus der 
Matro:;^"' ^P^'-^— Speichern der Meh.ahl von Messdaten ^n e^nem k^p^e rtt 

des.komprimie.en Mat_ durch.hrt. urn 

Vorrichtung nach Anspruch 8. wobei das Volumen eine H6he in dem System hat. wobei die statistlsche Analvse 
"men 2 ^-^^^''^'^'^e-achungsinfom.ation das zu dem S^ geh^^ge t 



Vorrichtung nach Anspruch 8, wobei 
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a) das System eine Vielzahi von Tanks aufweist, Oder 

b) das System einen untierirdischen Speichertank aufweist, oder 

c) das System einen oberirdischen Speichertank aufweist, oder 

d) das System einen teilweise oberirdischen Speichertank aufweist, oder 

e) die Messdaten gleichzeitig von den verschiedenen l\/lessvorrichtungen gesammelt werden, oder 

f) die von den Messvorrichtungen erhattenen Messdaten zu dem Computer ubertragen werden, oder 

g) die Vorrichtung weiter einen Host-Prozessor, der entfernt von dem Computer ist, aufweist, um die statisti- 
sche Analyse durchzufuhren, 

11. Vorrichtung nach Anspruch 8, wobei die Messvorrichtuhg ein Voiumenmessinstrument, eine Ausgabevorrichtung 
und eine Verkaufsaufzeichnungseinrichtung aufweist und vorzugsweise 

a) die Ausgabevorrichtung einen Totaiisierer umfasst oder 

b) die Verkaufsaufzeichnungselnrichtung den Betrieb eines Kassenterminals simuliert, oder 

c) das Voiumenmessinstrument einem in dem System angeordneten Sensor aufweist und der Sensor vor- 
zugsweise eine magnetostriktiver Tanksensor ist. 

12. Vorrichtung nach Anspruch 8, die weiter einen in dem System in Kontakt mit dem Volumen angeordneten Tempe- 
ratursensor aufweist, um eine Temperaturmessung des Volumens zu erhalten, und wobei vorzugsweise die Tem- 
peratumriessung dazu venwendet wird, um das Volumen zu berechnen. 

13. Vorrichtung nach Anspruch 8, wobei 

a) die Vorrichtung weiter eine Mehrzahl von an verschiedenen Orten innerhalb des Systems angeordneten 
Temperatursensoren aufweist, wobei jederTemperatursensor in Kontakt mit dem Volumen ist, um eine Mehr- 
zahl von Temperaturmessungen des Volumens zu erhalten, und wobei die Temperaturmessungen vorzugs- 
weise dazu verwendet werden, das Volumen zu berechnen, oder 

b) das komprimierte Matrixfomnat als Produkt einer Datenmatrix und der Transponierten der Datenmatrix er- 
zeugt wird, wobei vorzugsweise das Produkt durch Addition von Tellprodukten jeder aus einer Mehrzahl von 
Teiimatrizen der Datenmatrix mit der Transponierten der Teiimatrix gebildet wird. 

Revendlcatlons 

1. Proc6d6 de surveillance d'un systfeme de stockage et de distribution de fluide, ledit syst^me comprenant un ap- 
pareil de mesure pour mesurer un volume associe audit systeme, ledit precede comprenant : 

la collecte d'une plurality de donnees de mesure provenant dudit apparell de mesure sous une fornie exploi- 
table par un ordlnateur. ; 

le stockage de ladite pluralite de donnees de mesure sous un format matriciel compress6 dans une memoire 
d'ordinateur ; et 

I'analyse statistique dudit fonnat matriciei compress^ pour determiner une infomnation relative ^ la surveillance 
du fonctionnement. 

2. Precede seion la revendication 1 comprenant en outre les operations consistant^ 

a) determiner la presence de d^fauts de tonctionnementdans ledit systeme ; ou 

b) surveiller la precision dudit appareil de mesure ; ou . 

c) determiner le volume de fluide dans ledit systfeme ; ou 
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d) d6termlner si le retrait d'une certaine quantlt6 de flulde dudit systfeme est caus4 oar „n« f .it ^ ■ 
t.ma et de pr^fyence. d.livrer un signal avertissant qu'^ne fuite aT/d^I^^^^^^^^^^ IV'" 

e) rnterroger ledit appareil de mesure sous la commande dudit ordinateur • ou ^ ' 

svTtImT.'JHirTl^'".'^ • '^'^'^ ^ P««'^ temperature dudit volume dans ledit 

systems, led. caicul incluant de pr6f6rence la determination d'une valeur de correction a oatir Tl mZlt^l 
pond6r6e de la temperature dudit volume mesure en una pluralite d'endroltsTS^eu^Sit .J^^^^^^ 

sS^e'"'' ' '-'^te etape d-analyse 

3. Proc6d§ selon la revendlcation 1 , dans lequel 

c) ledit systeme comprend une pluralite de reservoirs ; ou ' ' 

d) ledit systfeme comprend un reservoir de stocl<age souten-ain • ou 

e) ledit systeme comprend un reservoir de stockage aerien • ou 

f) ledit systfeme comprend un reservoir de stockage partlelle'ment aerien • ou 

9) ladite etape de stockage comprend I'obtention dudit format matriciel compressfe en tant aue orodult d-..n« 

t^rlmlm^sTelti"^^^^ '"'"^ une jauge voiumfetrique, un appareil distributeur et un dispos.if d^enregis- 
i) ladite etape de coilecte comprend I'evaluation d'un volume initial de flulde dans ledit systfeme. 

c) evaiuer une valeur initiate dudit volume pendant ladite etape d'analyse 

b) ladte 6tape de oolteMs n'M pas «alI,4e lor.qu8 tediniulde Mt aJouM audits^ 
parei. de mesure destine . mesurer ladite haute'ur pour chacun'disdits vo[^^^^^^^^^^^ """^""""^ 



EP 0 904 526 B1 

Appareil pour surveiller un systeme de stockage et de distribution de fiuide, ledit appareil comprenant : 
un appareil de mesure destine k mesurer un volume associ6 audit systeme ; et 

un ordinateur comprenant un moyen de traltement pour collecter une pluralrte de donn^es de mesure aupres 
dudit appareil de mesure et une m6moire pour stocker (adite pluralite de donnees de hnesure sous un fonnat 
matriciel compress^ ; 

dans lequel ledit moyen de traitement r6alisant une analyse statfstique dudit fonnat matriciel compress6 en vue 
de determiner une infonnation relative a la surveillance du fonctionnement. 

Appareil selon la revendication 8, dans lequel ledit volume a une certaine hauteur dans ledit systeme, dans lequel 
ladite analyse statistique est une analyse par regression et dans lequel ladite information relative ^ la surveillance 
du fonctionnement est ledit volume associ6 audit systeme. 

Appareil selon la revendication 8, dans lequel 

a) ledit systeme comprend une pluralit6 de reservoirs ; ou 

b) ledit systeme comprend un reservoir de stockage souterrain ; ou 

c) ledit systeme comprend un reservoir de stockage aerien ; ou 

20 d) ledit systeme comprend un reservoir de stockage partiellement aerien ; ou 

e) ladite donnee de mesure est collectee simultanement aupres desdits divers appareils de mesure ; ou 

f) lesdites donnees de mesure obtenues aupres desdits, appareils de mesure sont transmises audit ordinateur ; 
ou 

g) ledit appareil comprend en outre un calculateur central se trouvant a distance dudit ordinateur pour realiser 
25 ladite analyse statistique. 

1 1 . Appareil selon la revendication 8. dans lequel ledit appareil de mesure inclut une jauge volumetrique, un appareil 
distributeur et un dispositif d'enregistrement des ventes et de preference 

30 a) ledit appareil distributeur comprend un totaiisateur ; ou 

b) ledit dispositif d'enregistrement des ventes simule le fonctionnement d'un temninal de lieu.de vente ; ou 

c) ladite jauge volumetrique comprend une sonde agencee dans ledit systeme et, de preference, ladite sonde 
est une sonde de reservoir a magnetostriction. 

35 12. Appareil selon la revendication 8 comprenant en outre un capteur de temperature, agence t I'interieur dudit sys- 
teme, en contact avec iedit volume pour obtenir une mesure de la. temperature dudit volume, et de prefefence 
dans lequel ladite mesure de temperature est utilisee pour calculer ledit volume, 

13. Appareil selon la revendication 8, dans lequel 

40 

a) I'appareil comprend en outre une pluralite de capteurs de temperature agences a differents endroits e 
I'interieur dudit systeme, chacun desdits capteurs etant en contact avec ledit volume pour obtenir une pluralite 
de mesures de temperature dudit volume, et de preference dans lequel lesdites mesures de temperature sont 
utilisees pour calculer ledit volume ; ou 
45 b) ledit format matriciel compresse est obtenu en tant que produit d'une matrice de donnees et de la transposee 

de ladite matrice de donnees et, de preference, dans lequel ledit produit est fomne par addition desdits produits 
partiels de chacune d'une pluralite de partitions de ladite matrice de donnees avec la transposee de chaque 
dite partition. 
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